A polygenic basis for birth weight in a wild population of red deer (Cervus elaphus).

The genetic architecture of traits under selection has important consequences for the response to selection and potentially for population viability. Early QTL mapping studies in wild populations have reported loci with large effect on trait variation. However, these results are contradicted by more recent genome-wide association analyses, which strongly support the idea that most quantitative traits have a polygenic basis. This study aims to re-evaluate the genetic architecture of a key morphological trait, birth weight, in a wild population of red deer (Cervus elaphus), using genomic approaches. A previous study using 93 microsatellite and allozyme markers and linkage mapping on a kindred of 364 deer detected a pronounced QTL on chromosome 21 explaining 29% of the variance in birth weight, suggesting that this trait is partly controlled by genes with large effects. Here, we used data for more than 2,300 calves genotyped at >39,000 SNP markers and two approaches to characterise the genetic architecture of birth weight. First, we performed a genome-wide association (GWA) analysis, using a genomic relatedness matrix to account for population structure. We found no SNPs significantly associated with birth weight. Second, we used genomic prediction to estimate the proportion of variance explained by each SNP and chromosome. This analysis confirmed that most genetic variance in birth weight was explained by loci with very small effect sizes. Third, we found that the proportion of variance explained by each chromosome was slightly positively correlated with its size. These three findings highlight a highly polygenic architecture for birth weight, which contradicts the previous QTL study. These results are probably explained by the differences in how associations are modelled between QTL mapping and GWA. Our study suggests that models of polygenic adaptation are the most appropriate to study the evolutionary trajectory of this trait.

Maternal effects do not resolve the paradox of stasis in birth weight in a wild red deer populaton.

In natural populations, quantitative traits seldom show short-term evolution at the rate predicted by evolutionary models. Resolving this "paradox of stasis" is a key goal in evolutionary biology, as it directly challenges our capacity to predict evolutionary change. One particularly promising hypothesis to explain the lack of evolutionary responses in a key offspring trait, body weight, is that positive selection on juveniles is counterbalanced by selection against maternal investment in offspring growth, given that reproduction is costly for the mothers. Here, we used data from one of the longest individual-based studies of a wild mammal population to test this hypothesis. We first showed that despite positive directional selection on birth weight, and heritable variation for this trait, no genetic change has been observed for birth weight over the past 47 years in the study population. Contrarily to our expectation, we found no evidence of selection against maternal investment in birth weight-if anything, selection favors mothers that produce large calves. Accordingly, we show that genetic change in birth weight over the study period is actually lower than that predicted from models including selection on maternal performance; ultimately our analysis here only deepens rather than resolves the paradox of stasis.

The role of maternally transferred antibodies in maternal performance in red deer.

Maternal effects are ubiquitous. Yet, the pathways through which maternal effects occur in wild mammals remain largely unknown. We hypothesise that maternal immune transfer is a key mechanism by which mothers can affect their offspring fitness, and that individual variation in maternally derived antibodies mainly depends on a mother's characteristics and the environmental conditions she experiences. To test this, we assayed six colostrum-derived antibodies in the plasma of 1447 neonates in a wild red deer population. Neonatal antibody levels were mainly affected by maternal genes, environmental variation and costs of prior reproductive investment. We found consistent heterogeneity in maternal performance across traits, with mothers producing the heaviest calves also having calves with more antibodies. Unexpectedly, antibody levels were not associated with calf survival. We provide a unique example of how evolutionary theory on maternal effects can be used to gain insight into the causes of maternal effects in wild populations.

The genetic architecture of maternal effects across ontogeny in the red deer.

Maternal effects, either environmental or genetic in origin, are an underappreciated source of phenotypic variance in natural populations. Maternal genetic effects have the potential to constrain or enhance the evolution of offspring traits depending on their magnitude and their genetic correlation with direct genetic effects. We estimated the maternal effect variance and its genetic component for 12 traits expressed over the life history in a pedigreed population of wild red deer (morphology, survival/longevity, breeding success). We only found support for maternal genetic effect variance in the two neonatal morphological traits: birth weight ( hMg2 = 0.31) and birth leg length ( hMg2 = 0.17). For these two traits, the genetic correlation between maternal and direct additive effects was not significantly different from zero, indicating no constraint to evolution from genetic architecture. In contrast, variance in maternal genetic effects enhanced the additive genetic variance available to respond to natural selection. Maternal effect variance was negligible for late-life traits. We found no evidence for sex differences in either the direct or maternal genetic architecture of offspring traits. Our results suggest that maternal genetic effect variance declines over the lifetime, but also that this additional heritable genetic variation may facilitate evolutionary responses of early-life traits.

Where is the optimum? Predicting the variation of selection along climatic gradients and the adaptive value of plasticity. A case study on tree phenology.

Many theoretical models predict when genetic evolution and phenotypic plasticity allow adaptation to changing environmental conditions. These models generally assume stabilizing selection around some optimal phenotype. We however often ignore how optimal phenotypes change with the environment, which limit our understanding of the adaptive value of phenotypic plasticity. Here, we propose an approach based on our knowledge of the causal relationships between climate, adaptive traits, and fitness to further these questions. This approach relies on a sensitivity analysis of the process-based model phenofit, which mathematically formalizes these causal relationships, to predict fitness landscapes and optimal budburst dates along elevation gradients in three major European tree species. Variation in the overall shape of the fitness landscape and resulting directional selection gradients were found to be mainly driven by temperature variation. The optimal budburst date was delayed with elevation, while the range of dates allowing high fitness narrowed and the maximal fitness at the optimum decreased. We also found that the plasticity of the budburst date should allow tracking the spatial variation in the optimal date, but with variable mismatch depending on the species, ranging from negligible mismatch in fir, moderate in beech, to large in oak. Phenotypic plasticity would therefore be more adaptive in fir and beech than in oak. In all species, we predicted stronger directional selection for earlier budburst date at higher elevation. The weak selection on budburst date in fir should result in the evolution of negligible genetic divergence, while beech and oak would evolve counter-gradient variation, where genetic and environmental effects are in opposite directions. Our study suggests that theoretical models should consider how whole fitness landscapes change with the environment. The approach introduced here has the potential to be developed for other traits and species to explore how populations will adapt to climate change.

Microgeographic adaptation and the effect of pollen flow on the adaptive potential of a temperate tree species.

In species with long-distance dispersal capacities and inhabiting a large ecological niche, local selection and gene flow are expected to be major evolutionary forces affecting the genetic adaptation of natural populations. Yet, in species such as trees, evidence of microgeographic adaptation and the quantitative assessment of the impact of gene flow on adaptive genetic variation are still limited. Here, we used extensive genetic and phenotypic data from European beech seedlings collected along an elevation gradient, and grown in a common garden, to study the signature of selection on the divergence of eleven potentially adaptive traits, and to assess the role of gene flow in resupplying adaptive genetic variation. We found a significant signal of adaptive differentiation among plots separated by < 1 km, with selection acting on growth and phenological traits. Consistent with theoretical expectations, our results suggest that pollen dispersal contributes to increase genetic diversity for these locally differentiated traits. Our results thus highlight that local selection is an important evolutionary force in natural tree populations and suggest that management interventions to facilitate movement of gametes along short ecological gradients would boost genetic diversity of individual tree populations, and enhance their adaptive potential to rapidly changing environments.

Tree, sex and size: Ecological determinants of male vs. female fecundity in three Fagus sylvatica stands.

Interindividual variation in fecundities has major consequences on population evolutionary potential, through genetic drift and selection. Using two spatially explicit mating models that analyse the genotypes of seeds and seedlings, we investigated the variation of male and female fecundities within and among three European beech (Fagus sylvatica) stands situated along an altitudinal gradient. Female and male individual fecundity distributions were both skewed in this monoecious species, and we found a higher variance in female as compared to male fecundities. Both female and male fecundities increased with tree size and decreased with density and competition in the neighbourhood, the details of these effects suggesting sex-specific strategies to deal with the impact of limited resource on fecundity. The studied populations were functionally male-biased. Among-individual variations in functional gender were not driven by tree size but by density and competition in the neighbourhood. Femaleness decreased under limited resource availability, an expected consequence of the higher cost of female reproduction. Considering the variation of gene flow and genetic drift across elevation, our results suggest that the adaptive potential could be enhanced by low genetic drift at low elevation, and by high pollen-mediated gene flow at high elevation. Finally, this study predicts a more efficient response to selection for traits related to male vs. female fitness, for a given selection intensity.

Ecological determinants of mating system within and between three Fagus sylvatica populations along an elevational gradient.

Studies addressing the variation of mating system between plant populations rarely account for the variability of these parameters between individuals within populations, although this variability is often non-negligible. Here, we propose a new direct method based on paternity analyses (Mixed Effect Mating Model) to estimate individual migration (mi ) and selfing rates (si ) together with the pollen dispersal kernel. Using this method and the KINDIST approach, we investigated the variation of mating system parameters within and between three populations of Fagus sylvatica along an elevational gradient. Among the mother trees, si varied from 0% to 48%, mi varied from 12% to 86% and the effective number of pollen donors (Nepi ) varied from 2 to 364. The mating patterns differed along the gradient, the top population showing higher m and lower s, and a trend to higher Nep than the bottom populations. The phenological lag shaped long-distance pollen flow both within population (by increasing mi at mother-tree level) and between populations (by increasing m at high elevation). Rather than the mate density, the canopy density was detected as a major mating system determinant within population; it acted as a barrier to pollen flow, decreasing the proportion of long-distance pollen flow and increasing si . Overall, the effects of ecological factors on mating system were not the same within vs. between populations across the gradient, and these factors also differed from those traditionally found to shape variation at range-wide scale, highlighting the interest of multiscale approaches.